Multiple hydrodynamical shocks induced by Raman effect in photonic crystal fibres 
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We theoretically predict the occurrence of multiple hydrodynamical-like shock phenomena in the 
propagation of ultrashort intense pulses in a suitably engineered photonic crystal fiber. The shocks 
are due to the Raman effect, which acts as a nonlocal term favoring their generation in the focusing 
regime. R is shown that the problem is mapped to shock formation in the presence of a slope and 
a gravity-like potential. The signature of multiple shocks in XFROG signals is unveiled. 
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I. INTRODUCTION 

The existence of regimes where nonlinear optical prop- 
agation and Bose-Einstein condensation (EEC) obey 
hydrodynamical-like models is well knowr>ir— . Among 
the resulting hydrodynamic-like phenomena for intense 
Hght pulses or spatial beams, the most fascinating is cer- 
tainly the formation of shocks^, which can also be con- 
sidered in Bosc-Einstein condensates (BEC)^"^^. 

Theoretical and experimental work on optical shock 
formation has been recently reported for nonlocal me- 
dia, like thermal liquidsi^""— , liquid crystals^'^ and pho- 
torefractive material s ^^'^° . In the time domain, shocks 
have been considered in connection with the carrier 
wave^i and quadratic media3^. Spatiotemporal hyper- 
bolic shocks have also been studied^^, together with re- 
lated nonlinear X-wave generatiort2ii2^ and multidimen- 
sional effects^^iSl. 

When considering temporal shocks, if it is true that 
the development of extreme nonlinear optics^^ unavoid- 
ably requires the consideration of ultra-wide spectral- 
band and intense excitations (common to all shock phe- 
nomena), the presence of higher order effects radically 
affects shock formation and related regularization pro- 
cesses, such as undular bore a^^'^° . In this respect, mi- 
crostructured photonic-crystal fibers (PCFs)^ offer an 
unprecedented framework for exploring the hydrodynam- 
ical properties of light. Indeed, not only do solid-core 
PCFs exhibit very pronounced nonlinear effects (mainly 
thanks to their small effective modal areas), but their 
dispersion profile can be engineered to a large degree. 
The latter circumstance is particularly appealing for the 
physics of shock wave generation, because most of the dy- 
namics described by nonlinear Schrodinger (NLS) models 
predict wave-breaking phenomena in the limit of vanish- 
ing dispersion. Such a condition, for temporal shocks, 
is in general very difficult to achieve over wide spectral 
regions. PCFs can however be designed with almost flat 
dispersion over bandwidths of several hundreds nanome- 
ters, thus opening up new perspectives in temporal opti- 
cal shock generation and related applications. 

From a theoretical point of view, there are still several 
open problems concerning the physics of optical shocks. 
One of the most challenging is the onset of wave-breaking 
phenomena in the regime of focusing nonlinearity and 



anomalous group velocity dispersion. Even if unstable 
kink-antikink solutions (regularized shock fronts with- 
out oscillations) are known to exist when including Ra- 
man term&'^^'?^=, the standard hydrodynamical limit of the 
NLS apparently rules out shock formation in the focus- 
ing regime, as the jump condition on the resulting Euler 
equations (see, e.g,!!^) gives an imaginary result, which 
does not have a simple physical interpretation. However, 
it has been shown that, in non-local modelsi^ for spatial 
beam propagation and BEC, a multi- valued solution of 
the hydrodynamical equation for the so-called velocity 
held can be predicted by using the method of charac- 
teristics. Such a solution corresponds in the numerical 
simulations of the nonlocal NLS to the clear formation 
of an undular-like bore in a regime in which it competes 
with the modulational instability, as also considered in—. 
The scenario is hence extremely rich and interesting, and 
the open question is whether temporal shocks, undular 
bores and wave-breaking phenomena may form in a fo- 
cusing regime when considering a real-world system. 

The most natural counterpart of a nonlocal nonlinear 
response in the time domain is the Raman effect. Indeed 
this is described by a kernel response function, which, un- 
der suitable conditions on the pulse duration, leads to a 
nonlinear shock term containing the first-order derivative 
of the intensitySi. Previous investigations of intense light 
propagation in solid-core PCFs when including the Ra- 
man effect have shown that the peculiar breathing phe- 
nomenon exhibited by higher-order solitons can be used 
to excite the formation of quasi-symmetric resonant ra- 
diation in a step-like fashion in presence of highly dis- 
torted GVD^. However, in that work, and in all pre- 
vious observations of the phenomenon (see e.g;^), there 
is no convincing explanation of why the soliton breath- 
ing should increase its rate in the presence of the Raman 
effect, which leads to soliton splitting, and the internal 
dynamics of the breathing has currently a rather cum- 
bersome interpretation. 

In this paper we show that under suitable conditions 
the soliton breathing process is accompanied by the for- 
mation of multiple optical shock waves leading to strong 
spectral broadening, and that this can described by a hy- 
drodynamical model containing a gravity-like slope term, 
or equivalently a constant external electric field in a cold 
plasma. This point of view allows one to clarify and 
shed new light on the several well-known phenomena de- 



scribed above. Our theoretical results are validated by 
real-world simulations in a realistic PCF, and represent 
the first theoretical prediction of shock waves in the fo- 
cusing regime, also unveiling their measurable signatures 
in the XFROG signal. 

The paper is organized as follows. In section |TT] we 
discuss the hydrodynamical model for the NLS in the 
presence of a Raman term. In section IIIII we adopt the 
method of characteristics to predict and give a mechani- 
cal interpretation to the occurrence of multi- valued solu- 
tions for the velocity- field. In section HVl we describe the 
PCF we used for investigating temporal hydrodynamic- 
like shocks and report on the corresponding simulations 
and the XFROG signal at the shock formation. Conclu- 
sions are drawn in section |Vl 



II. THE HYDRODYNAMICAL LIMIT IN THE 
ACCELERATED REFERENCE FRAME 

We start by considering the simplest model for an op- 
tical pulse propagating in an optical fibre, described by 
the envelope equation^! 
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In Eq. ([TJ, ijj is the envelope of the electric field, scaled 
with the soliton power Pq — |/32|/(7io), where (32 is the 
second order dispersion coefficient, 7 is the nonlinear co- 
efficient of the fiber, and to is the input pulse duration. 
The dimensionless propagation length z is scaled with the 
second order dispersion length Ld2 — ^o/ 1/^2 1, while the 
temporal variable t is scaled with to. The last term in ([T]) 
represents the Raman effect, and tj^ = Tj^/to, where Tj^ 
is the Raman response time (about 2 fs in silica)^'*. For 
simplicity, only second order dispersion has been included 
in the model of Ec^. ([1]). Note that Eq. (HI) is written in 
the anomalous dispersion regime, where bright solitons 
are expected in the absence of the Raman term (th — 0). 
Eq. ([1]) is subject to the initial condition 



^(z = 0,t) =iVsech(i). 



(2) 



The generation of the shock-like dynamics corresponds 
to the case of large N, and introducing the smallness pa- 
rameter e ~ I/VN, the hydrodynamical limit is obtained 
via rewriting Eq. ([T]) in the semiclassical form with the 
scaled variables z — > z/e^, t — > t/e^, ip — >■ ipe, obtaining 

2 

JeV'. + yV'« + |V'|'V'-^MI^I' = 0. (3) 

The Raman term appears at the relevant order in the 
hydrodynamical limit for r ~ e, which is the condition 
at which the Raman effect radically alters the shock dy- 
namics, as also confirmed by the simulations reported 
below. 

It is interesting to derive the hydrodynamical limit in 
the accelerated reference frame where the soliton moves 



in presence of the Raman effect^. This is obtained by 
introducing the variable ^ = t—^z^ (with g = 32r/{a^/15 
as in^), and performing the Gagnon-Belanger phase 
transformation^^ 

„2 

i> = f{tz)exp[-i{^z^^gzt + az)], (4) 

after which the resulting NLS is written as^i^ 

^fz - gU + af + i/^^ + I/IV - TRfd^lff = (5) 

Letting / — Y^exp(i0), and performing a WKB expan- 
sion on Eq. dS]), one obtains 



p,+d^{pv)^0 (6) 



i^| = -.^ + a + p-r.a,p+i-L|^ (7) 



By deriving Eq. ([7]) with respect to £,, one has (after 
defining the velocity field v = (p^, physically correspond- 
ing to the instantaneous frequency inside the pulse) 

v, + vv^^-d^{UQP + U). (8) 

The quantum pressure potential is defined as 
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and the modified 'classical' potential term is 

U^g^-p + TRd^p. (10) 

Note that potential U in Eq. pH)) is exactly the same 
as the one formerly introduced in"^^. The above relations 
show that the hydrodynamical limit of the NLS with Ra- 
man corresponds to fluid motion in an accelerated frame 
(being g the acceleration) or, equivalently, to a charged 
plasma in a constant electric field of magnitude g. 



HI. EFFECTIVE PARTICLE 
INTERPRETATION 

Eq. ([8]) can be solved before the occurrence of shocks 
in an approximate way by recalling that the equation 
for p is obtained at a higher order in e, and hence that 
the initial dynamics is ruled by the velocity field v only. 
This allows, as a zero-order approximation, to take p(^) 
as a fixed profile, e.g. as a sech function denoting the 
initial pulse profile. Eq. ([S]) can hence be mapped by 
the method of characteristics into the system of ordinary 
differential equations 
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dz 
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(11) 



with the initial conditions ^ = s and v = a.t z — 
0. The solution for a given /o(^), gives the manifold 



(^(s, z),v{s, z)), which for a given z provides a parametric 
pfot of V versus z and aUows one to predict the occurrence 
of hydrodynamical shocks. 

Eqs. (|lip can be interpreted as the motion of a par- 
ticle with trajectory ^(z) in the potential U and can be 
rewritten as a single Newton-like equation 



d^ 



dU 



(12) 



Shocks correspond to the case in which the plot of the 
trajectory in the space (^, v) display a vertical slope, i.e. 
d^/dv = 0. 

This can be geometrically interpreted as follows; var- 
ious trajectories are generated by varying the particle's 
initial position ^ = s, see Fig. [TJa), which corresponds 
to considering the motion of particles in the potential U 
[shown in Fig. [ijb)] for different values of tr which are 
falling (with initially zero velocity) from position ^ = s. 
A shock, i.e. a multi-valued function v = w(^) in the 
plane {£,,v) [see Fig. HJc)] corresponds to the existence 
of trajectories reaching the same position ^ at the same 
z with different velocities, i.e. to collisions between parti- 
cles falling from different initial positions under the effect 
oiU. 

Given the shape of U{£,), such collisions specifically 
involve the particles trapped in a bounded motion, which 
implies the formation of shocks in a periodic fashion, see 
also Fig. [TJa) and (c). 

As Tn increases the bounded motion occurs within an 
increasingly smaller region in ^ [the potential well is re- 
duced, as it is clear from Fig. [Ifb)], and it becomes 
negligible for larger values of tr. This implies that for 
increasing r^j, shocks become more frequent (as the cor- 
responding particle collisions), but also that the corre- 
sponding discontinuities will be less pronounced. Such a 
qualitative interpretation, gives a simple picture of what 
is observed in the simulations discussed below. 



IV. SIMULATIONS 

In this section we show that the hydrodynamical limit 
of the NLS equation [see Eq. dSj] can be realized in 
practice through a realistic PCF design. We also compare 
the simulation dynamics of the idealized model of Eq. 
([3]) with a more realistic simulation based on the full 
generalized NLS equation (GNLSE). 

We start by integrating numerically Eq. ([3]) with pa- 
rameters Tn — 0.1 and e = 0.2, with initial condition 
Tpin = sech (t) . This set of parameters is chosen to re- 
alize the hydrodynamical limit by means of the scaling 
discussed in section HIl After a propagation length of 
z ~ 0.62, one can observe the formation of the first shock, 
see Fig. [2Ja) where the time-derivative of the phase (i.e. 
the velocity field) v{t) = cpt has been plotted. Linked 
to the formation of the phase shock, one observes the 
formation of a sharp cusp in the time domain located 
at the maximum amplitude of the pulse, see Fig. [Hb). 
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FIG. 1: (Color online) (a) Positions of a particle falling inside 
potential U{^) as described by Eq. (fTO|l . for different values 
of the initial position ^ = s and for tr = 0.1. Multiple col- 
lisions between oscillating bound particles produce multiple 
values for the velocity field v, i.e. shocks. Non-bound par- 
ticles can also collide with bound particles - although more 
rarely - contributing to isolated shocks, (b) Shape of the ef- 
fective potential for different values of tr with p = sech (^) ; 
(c) Solution of the hydrodynamical system (jlip for different 
values of z, showing multiple-shock occurrences. The velocity 
is plotted for clarity with an arbitrary vertical shift for the 



various z. 



Large values of the time-derivative at the cusp produce 
large spectral broadening in the frequency domain, which 
in the conventional and very well-known picture is asso- 
ciated with the temporal and spectral 'breathing' oscil- 
lations of higher-order solitons^li^^. In the presence of 
the Raman effect (t^ 7^ 0), one will observe the occur- 
rence of multiple shocks, which become more and more 
frequent during propagation, contrary to the perfectly 
regular oscillations of the soliton breathing phenomenon 
in the absence of the RamanSiiS. This 'breathing accel- 
eration' can again be interpreted as multiple shocks of the 
falling fluid induced by the 'gravity-like' potential of Eq. 
([5]), as described in section Hill Moreover, as is shown 
in Fig. mja), a very clear undular bore phenomenon de- 
velops in the trailing edge of the Raman-shifting pulse 
at longer propagation distances, which does not occur 
in the absence of the Raman effect. This is associated 
to the strong oscillations of the velocity field v induced 
by the quantum pressure potential ^. This term be- 
comes indeed important in proximity of those parts of 
the pulse for which p — >■ 0, i.e. far from the pulse cen- 
ter. This, when combined with the Raman gravity-like 
potential PH)) . leads to the formation and transport of 
increasingly stronger oscillations in the trailing edge of 
each pulse. 




FIG. 2: (Color online) Pulse propagation and shock forma- 
tion according to Eq. (I13[) . using the GVD of Fig. [S] (a) 
Velocity field v = dt<j) as a function of time for three different 
propagation distances, namely z = 0.4 (red solid line, before 
the occurrence of shock), z = 0.62 (black dotted line, at the 
moment of shock) and z = 1.5 (blue dashed line, long distance 
dynamics), for tr — 0.1 and e — 0.2. For long propagation 
dynamics, the formation of clear undular bores is observed, 
(b) Amplitude \tp\ at the same values of z as in (a). At the 
moment of the first shock, a sharp cusp is formed. 




> -200 
CD 



-400 



600 800 1000 1200 
wavelength, nm 



1400 



FIG. 3: (Color online) GVD curve of the PCF used in our sim- 
ulations. The black dots represent the data measured exper- 
imentally through white-light interferometry, while the solid 
curve is a polynomial fit. The ZDW is close to Xz ~ 690 nm. 
The inset shows an SEM of the fiber cross section. 



We now compare the shock dynamics of the above sim- 
plified model with the full GNLSE, applied to a realistic 
waveguide. In Fig. [3] we show the GVD of a solid-core 
PCF, the SEM picture of which is given in the inset. 
The black dots correspond to measured data, while the 
red solid line shows a fit. The fiber has a core diameter 
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FIG. 4: (Color online) XFROG spectrograms for the pulses 
at z equal to (a) 1.761, (b) 2.192, (c) 2.723, and (d) 3.393 cm. 
The peak power of the to = 130 fs input pulse was P = 5 kW 
{N = 3.1), and they were launched at A = 1 /xm. Fig (a) and 
(b) are snapshots of the first shock. Fig (c) and (d) show the 
second shock, coinciding with the second spectral oscillation. 



of approximately 1.5 /im, an estimated nonlinear coeffi- 
cient 7 ~ 100 W~^m~^, and a zero-GVD point located 
at Xz ~ 690 nm. We use the following GNLSE^: 



id.i> + D{idt)ip - 






^ 



r{t--t')\^{t')\'dt' 



0, 



(13) 

where r{t) = [(1 — 0)S{t) + 9h{t)]/to is the total response 
function, to is the input pulse duration, 9 is the relative 
contribution between the non-instantaneous Raman and 
the instantaneous Kerr effect in the material (6* = 0.18 
in silica), S{t) is the Dirac-delta function. The Raman 
response function is introduced in the code by using the 
HoUenbeck and Cantrell approach^. The /i-function is 
normalized, /_ r{t')dt' = 1, while the dimensionless 
Raman response time parameter used in Eq. ([T|) and 
Eq. ([3]), proportional to the first momentum, is given by 
TR = {l/h)C^t'r{t')dt'. Operator b{idt) in Eq. ^ 
describes the full complexity of the fiber GVD shown 
in Fig. El and is given by b{idt) = Em>2 ^'"mt'' , 
where /?„ = {d"'' f3{u!) / duj™)uj=uio is the m-th derivative 
of the propagation constant /3{uj) calculated at a suitable 
arbitrary reference frequency ujq. 

The panels in Fig. |4|a-d) show the evolution of the 

XFROG trace 1(6, t) = \A{t)A,^f{t - T)e~'^*dtf^ (A.^t 
is a suitable reference pulse) recorded at four different 
propagation distances, corresponding to the first two oc- 
currences of maximum spectral broadening, for A = 1 
/im, relatively far from the zero GVD point. 

The goal of our simulations is to identify the multiple 
shock generation during pulse propagation in a PCF. In 
section II we showed that the hydrodynamical limit re- 
quires large N, thus very high peak powers. We find how- 



100 




FIG. 5: (Color online) (a-d) Amplitude (black dotted line) 
and velocity field (red solid line) in the same conditions as in 
Figs. [4ja-d). (a) is the moment of the first shock, while (b-d) 
show the generation of undular bores from the strong veloc- 
ity field oscillations of the singularities due to the quantum 
pressure potential ([9]). 



ever, that in our siumlations the dynamics of the multiple 
shock can best be demonstrated for lower N. In this case, 
the dynamics are not obscured by unwanted effects like 
resonant radiation, which disturb the soliton propagation 
and prevents thus a clear signature of the multiple shocks. 
In this regard, the presented GVD curve in Fig. [Hlwill be 
advanteous for the purpose of observing shock phenom- 
ena. A higher value of /?2 would lead to a massive Raman 
self-frequency shift which would obscure shock formation. 
On the other hand, smaller values of (52 (or too large 
values of N) lead to a dominating self-phase modula- 
tion, since Ld2 would be much longer that the nonlinear 
length Lnl = (T-P)"^- In that case, a very strong spec- 
tral broadening would be accompanied by the emission of 
resonant radiation from solitons"*^, which would also dis- 
turb the formation of clear shocks. The fiber is pumped 
by a sech pulse with N — 3.1. The amplitude shocks that 
develop at these moments are evident. The panels in Fig. 
[5]Ja-d) show the corresponding intensity profiles and the 
velocity field v{t) = dtcj) at the same values of z as in Fig. 
131 Again, one can see phase discontinuities located at 
the various pulse centers. Moreover, one can notice the 
generation of strong undular bores in the trailing edge of 
the pulse, which agree qualitatively with the simplified 
model examined above. 

In Figs. ini^a,b) we show the evolution of Eq. P^ in 
two different cases, one for t/j — 0.0334 [Fig. Hl^a)] and 
the other for th = 0.0557 [Fig. IH^b)]. One can notice 



the well-known slight increase of the rate of spectral os- 
cillations when increasing the value of t/j, which cannot 
find an easy explanation in the model of Eq. (|13p . How- 
ever, such behavior can easily be explained in terms of 
the hydrodynamical analogy expressed by Eqs. (I6][7l): in 
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FIG. 6: (Color online) Spectral evolution of P = 5 kW pulses 
for Raman fractions (a) fa = 0.18 (th = 0.0334), and (b) 
/h = 0.3 [tr — 0.0557), with to ~ 130 fs and a pulse launched 
at A = 800 nm in the PCF of Fig. (3] For increasing values of 
tr, the spectral oscillation rate rises, which can be interpreted 
by the hydrodynamical formulation of Eqs. (|6I7|) . Horizontal 
dashed lines indicate the distance at which one has 5 spectral 
oscillations, to help visualization. 



presence of Raman effect, the 'falling' photon fluid feels 
a gravity-like field which increases the rate of shocks dur- 
ing propagation, according to the qualitative explanation 
based on colliding particles given in section Hill 



V. CONCLUSIONS 

In conclusion we have shown that suitably engineered 
photonic crystal fibers may sustain the formation of mul- 
tiple hydrodynamical-like shocks during the propagation 
of ultrashort pulses. The shocks are shown to be clearly 
evident in the XFROG signal, have measurable signa- 
tures, and occur for a focusing nonlinearity. This effect 
can be exploited for the optimization and control of su- 
percontinuum generation, and is expected to play a role 
when considering multidimensional dynamics. 
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